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Some recent theoretical work on inclusive B decays relevant for the model independent determi- 
nation of \Vub\ and |Vcf,| is summarized. The theoretical predictions and their reliability for several 
differential decay distributions in B — * Xc,ueu and B — > Xs'y are reviewed. These can be used to 
determine certain important HQET matrix elements. The upsilon expansion and ways of testing it 
are discussed. 



I. INTRODUCTION 

In the near future, a large part of the high energy experimental program will be devoted to testing the Cabibbo- 
Kobayashi-Maskawa (CKM) picture of quark mixing and CP violation by directly measuring the sides and (some) 
angles of the unitarity triangle. If the value of sin(2/3), the CP asymmetry in S ^ J/ipKs, is not too far from the 
CDF central value [|l| , then searching for new physics at the B factories will require a combination of several precision 
measurements. Particularly important are \Vub\ and |Vtd|, which are the least precisely known elements of the CKM 
matrix. The latter will be measured hopefully in the upcoming run of the Tevatron from the ratio of Bg and B^ 
mixing, and will not be discussed here. This talk is motivated by trying to understand what the chances are to 

• ...reduce (conservative) error in \Vcb\ below ^ 5%? Although sometimes a smaller error is quoted already 
(e.g., by the Particle Data Group), it is hard to bound model independently the possible quark-hadron duality 
violation in the inclusive, and the size of ^ corrections in the exclusive determination. 

• ...determine \Vub\ with less than 10% error? The inclusive measurements require significant cuts on the 
available phase space; the exclusive measurements require knowledge of the form factors. 

«... reduce the theoretical uncertainties in the B — > Xg^ photon spectrum? The effect of the experimental cut 
on the photon energy needs to be better understood, and there are subtleties in the OPE beyond leading order. 
The theoretical reliablility of inclusive measurements can be competitive with the exclusive ones (or even better in 
some cases). For example, for the determination of \Vcb\ model dependence enters at the same order of AQQj~,/m^ 
corrections from both the inclusive semileptonic B — s- XcCD width and the B — > D*ev rate near zero recoil. 

Inclusive B decay rates can be computed model independently in a series in Aqcd/^6, using an operator product 
expansion (OPE) |^-^. The rrib — s- oo limit is given by h quark decay. For most quantities of interest, this result is 
known including the order as and the dominant part of the order corrections. Observables which do not depend 
on the four- momentum of the hadronic final state (e.g., total decay rate and lepton spectra) receive no correction 
at order AgcD/'rif, when written in terms of rrih, whereas differential rates with respect to hadronic variables (e.g., 
hadronic energy and invariant mass spectra) also depend on K/nrib, where A is the vtlb — rub mass difference in the 
nib — *■ oo limit. At order AQ(-.j-,/m^, the corrections are parameterized by two hadronic matrix elements, usually 
denoted by Ai and A2. The value A2 — 0.12GeV^ is known from the B* ~ B mass splitting. 

For inclusive b q decay, corrections to the nib 00 limit are expected to be under control in parts of 
phase space where several hadronic final states are allowed to contribute with invariant masses satisfying ?7i^ > 
-I- (few times) AqcDm-fc-Q Such observables which "average" sufficiently over different hadronic final states may be 
predicted reliably. (We are just beginning to learn quantitatively what sufficient averaging is.) The major uncertainty 
in these predictions is from the values of the quark masses and Ai, or equivalently, the values of A and Ai. These 
quantities can be extracted from heavy meson decay spectra, which is the subject of a large part of this talk. 



*Invited talk at the American Physical Society, Division of Particles and Fields Conference, Jan. 5-9, 1999, Los Angeles, CA. 
'^However, it is not true that several such states are required to contribute: e.g., B — s- XcCU decay in the small velocity limit 
can be computed reliably, even though it is saturated by D and D* only Mpl. 
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An important theoretical subtlety is related to the fact that A cannot be defined unambiguously beyond perturbation 
theory ^ , and its value extracted from data using theoretical expressions valid to different orders in the as may vary 
by order Aqcd- However, these ambiguities cancel |^ when one relates consistently physical observables to one 
another, i.e., the formulae used to determine A and Ai contain the same orders in the perturbation series as those 
in which the so extracted values are applied to make predictions. 



II. THE HQET PARAMETERS A AND Ai 

The shape of the lepton energy [p[-pl| or hadronic invariant mass spectrum in semileptonic B — > Xcii^ 

decay, and the photon spectrum in i? — > p"^-p^ can be used to measure the heavy quark effective theory (HQET) 
parameters A and Ai. Testing our understanding of the B — s- XcCv spectra is important to assess the reliability of 
the inclusive determination of |Vc6|, and especially that of |Vu&|- Understanding the photon spectrum in _B ^ Xsj 
is important to evaluate how precisely the total rate can be predicted in the presence of an experimental cut on the 
photon energy . Studying these distributions is also useful to establish the limitations of models which were built 
to fit the lepton energy spectrum in semileptonic decay. 

The theoretical predictions are known to order a^Po, where po = 11 — 2n//3 is the coeflticient of the one- loop 
/3-function. This part of the order piece usually provides a reliable estimate of the full order correction, and it 
is straightforward to compute using the method of Smith and Voloshin |^ . In some cases the full correction is 
known. The order Aq^p/mj^ terms have also been studied ^1 13,1^], and are used to estimate the uncertainties. 



The OPE for semileptonic (or radiative) B decay does not reproduce the physical spectrum point-by-point in 
regions of phase space where the hadronic invariant mass of the final state is restricted. For example, the maximum 
electron energy for a particular hadronic final state X is Ee'^^'^^^ — (m^ — m\) /2mB, so comparison with experimental 
data near the endpoint can only be made after sufficient smearing, or after integrating over a large enough region. 
The minimal size of this region was estimated to be around 300 — 500 MeV Q . The hadron mass spectrum cannot 
be predicted point-by-point without additional assumptions, but moments of it are calculable. In general, higher 
moments of a distribution or moments over smaller regions are less reliable than lower moments or moments over 
larger regions. The strategy is to find observables which are sensitive to A and Ai, but the deviations from h quark 
decay are small, so that the contributions from operators with dimension greater than 5 are not too important. 



A. B XcCV decay spectra 

Last year the CLEO Collaboration measured the first two moments of the hadronic invariant mass-squared (sh) 
distribution, {sh — rnj^) and {{sh — m|,)^), subject to the constraint Ee > 1.5 GeV Here rno = {rriD -l-3m£)*)/4. 
Each of these measurements give an allowed band on the A — Ai plane. These bands are almost perpendicular, so 
they give the fairly small intersection region shown in Fig. 1. The central values at order are [ pl[ 

A = (0.33 ±0.08) GeV, Ai = -(0.13 ± 0.06) GeV^ (1) 

The unknown order ^^^^qj^ / rrvl corrections not included in this result introduce a large uncertainty, especially for 
the second moment. As a result, the allowed range is much longer in the direction of the first moment band than 
perpendicular to it; see Fig. 2. 

Similar information on A and Ai can be obtained from the lepton energy spectrum, dV/dEe- It has been measured 
both by demanding only one charged lepton tag, and using a double tagged data sample where the charge of a 
high momentum lepton determines whether the other lepton in the event comes directly from semileptonic B decay 
(primary) or from the semileptonic decay of a i? decay product charmed hadron (secondary). The single tagged 
data has smaller statistical errors, but it is significantly contaminated by secondary leptons below about 1.5 GeV. 
Therefore, Ref. M considered the observables 
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FIG. 1. Bands in the A — Ai plane defined by the mea- 
sured first and second moments of the hadronic mass- 
squared and lepton energy distributions. Note that 
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FIG. 2. Estimates of the theoretical uncertainty in A and Ai 
due to unknown l/m^ terms from the shape of the electron spec- 
trum (solid ellipse), from hadronic mass moments (long dashed 
ellipse), and the B — > Xs'y spectrum (short dashed ellipse; see 
Sec. IIB). Only the relative size and orientation of the ellipses 
are meaningful, not their position. (From Refs. jl^ and 
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I, which is in good 



Using the CLEO data the central values A = 0.39 GeV and Ai = -0.19 GeV^ were obtained _ 
agreement with Eq. (0). 

However, following Ref. [0, CLEO determined the first two moments of the spectrum, (Ee) and {{Ee — (i?e))^), 
without any restriction on Ee, using the double tagged data and an extrapolation to E^ < 0.6 GeV. The result of this 
analysis is also plotted in Fig. 1, and yields quite improbable values for A and Ai. The extrapolation to E^ < 0.6 GeV 
introduces unnecessary model dependence, so this result may be less reliable than the hadronic invariant mass analysis. 
The extracted values of A and Ai are very sensitive to small systematic effects; and it seems problematic for the 
exclusive models used for the extrapolation to simultaneously reproduce the inclusive lepton spectrum and the B 
semileptonic branching fraction p3t . 

There are a number of points to emphasize regarding these results: 

1. Taking A from Eq. gives a determination of \Vcb\ from semileptonic B width with ^ 3% uncertainty. 

2. The inclusive \ Vcb\ seems to be slightly larger than the exclusive (maybe not "significantly" but "consistently"). 
Can the theoretical uncertainties be reduced by combining the exclusive and inclusive determinations of iVcfcl? 

3. Since the bands from the lepton energy spectrum and the first moment of the hadronic mass-squared spectrum 
are almost parallel, an independent constraint on the A — Ai plane is needed. (See Section II B.) 

4. Models do not seem to do well for B D**ev and B —> Z?(*'7rez/ — relying on them may be dangerous. 
Composition of semileptonic B decay seems to be not really understood yet. (See next subsection.) 



1. Semileptonic B decays to excited charmed mesons 

In the heavy quark symmetry limit p^ , the spin and parity of the light degrees of freedom in a heavy meson are 
conserved. In the charm sector, the ground state is the (Z?, D*) doublet of heavy quark spin symmetry with spin- 
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parity of the light degrees of freedom sj" = . The four lightest excited states, sometimes referred to as D**, are the 

i"^ and the (Di, D^) doublet with sj* = |" 



{D^, Dl) doublet with sj' = and the {Di, D^) doublet with sj* = | . The Di and have been observed with 



masses near 2.42 and 2.46 GeV, respectively, and width around 20MeV. States in the s^' = ^ doublet can decay 
into D^*'^TT in an s-wave, and so they are expected to be much broader than the Di and Z?2 which can only decay in 
a d-wave. (An s-wave decay for the Di is forbidden by heavy quark spin symmetry ||25|| .) The first observation of the 
D\ state with mass and width about 2.46 GeV and 290 MeV, respectively, was reported at this Conference p6| . 

B Dieu and B D2ev account for sizable fractions of semileptonic B decays, and are probably the only 
three-body semileptonic B decays (other than B D''*^ev) whose differential decay distributions will be precisely 
measured. ALEPH and CLEO measured recently, with some assumptions, B{B — > DieD^) — (6.0 ± 1.1) x 10^"^ [ p7| , 
while B — > D2eDf, has not been seen yet. Heavy quark symmetry implies that in the mg — s- oo limit {Q = c, 6), matrix 
elements of the weak currents between a B meson and an excited charmed meson vanish at zero recoil. However, in 
some cases at order Aqcd/toq these matrix elements are not zero. Since most of the phase space for semileptonic B 
decay to excited charmed mesons is near zero recoil, 1 < v ■ v' ^ 1.3, Aqcd/^^q corrections can be very important. 

The matrix elements of the weak currents between B mesons and Di or D2 mesons are conventionally parameterized 
in terms of a set of eight form factors fi and ki ||2^ . At zero recoil only /y^ , defined by 

{D, («', e)\cr b \B{v)) = ^mo, tub [fv, e*^ + (/y.«'^ + Inv'") (e* • «)] , (3) 

can contribute to the matrix elements. In the mg — > 00 limit, fi and ki are given in terms of a single Isgur-Wise 
fmiction, t{w) Heavy quark symmetry docs not fix t(1), since fv^ — (1 — w^) t(w) +C'(1/toq), and so fvii^) = 
in the infinite mass limit independent of the value of t(1). 

At order l/mq several new Isgur-Wise functions occur, together with the parameter A', which is the analog of A 
for the (Di, D2) doublet. A' — A ~ 0.39 GeV follows from the measured meson masses |2^]. At order l/mg, fvi (1) is 
no longer zero, but it can be written in terms of A' — A and the Isgur-Wise function t{w) evaluated at zero recoil [ p8| , 

V6/,(l) = -i<^^.(l). (4) 

This relation means that at zero recoil heavy quark symmetry gives some model independent information about the 
l/mg corrections, similar to Luke's theorem |^o| for the decay into the ground state {D, D*) doublet. 

Since the allowed kinematic range for B Dieve and B Djei^e decay are fairly small {1 < w ^ 1.3), and there 
are some constraints on the 1 /mg corrections at zero recoil, it is useful to consider the decay rates expanded in powers 
of w — 1. The general structure of the expansion of dT/dw is elucidated schematically below, 



dw 



Vw^ — 1 



^ Vw^~l 

dw 

(|A|=04) 



(0 + Oe + + ...) + (w - 1) (0 + £ + ...) -I- (w - 1)2 (1 +£ + ...) + . . 
(0 + Oe + + ...) + (w - 1) (1 + £ + ...) + (w - 1)2 (1 + e + ...) + . . 



dw 



(u;2 _ ^)3/2 (1 + e + ...) + (u; _ i)(i + e ...) + ... . (5) 



Here A is the helicity of the Di or ZJj: and e" denotes a term of order (Aqcd/^tiq)"- The zeros in Eq. (||) are 
consequences of heavy quark symmetry. 

Table I summarizes some of the more important predictions. The first row (iJoo) shows the infinite mass limit. 
Approximations Bi and B2 are two different ways of treating unknown 1 /mg corrections, and may give an indication 
of the uncertainty at this order. One of the most interesting predictions is that the i? — > £>! decay rate should 
be larger than B —>■ D2, contrary to the infinite mass limit, but in agreement with the data. A number of other 
predictions for different helicity amplitudes, factorization of i? — )■ (I?i, D2)neD decay, sum rules, etc., are presented 
in Ref. ||2^ (see also |0). The main points are: 

1. At zero recoil, order l/mg contributions to semileptonic B —>■ Di, decays (any excited charmed meson with 
-|- parity) are determined by the mg — > 00 Isgur-Wise function and known hadron mass splittings. 
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Approximation 


R = To'/rDi 


r(l) 


6.0 X 10"^ 
B{B Diev^) 


1/2 


^ Di+D*+Dl + Dl /T Di 


Boo 


1.65 




1.24 




4.26 


Bi 


0.52 




0.71 




2.55 


B2 


0.67 




0.75 




2.71 



TABLE I. Predictions for Vd'/Td^, t(1), and V u^+di+di+d-^/^ d^- (From Ref. ||.) 



2. Decay spectra can be predicted near zero recoil, including the 1/toq corrections, with reasonable assumptions. 

3. Test heavy quark symmetry for B decays to excited charmed mesons, where the l/mg terms are sometimes the 
leading contributions. Large l/mg corrections to some predictions can be checked against data. 

4. Better understanding of inclusive — ^ exclusive in semileptonic B decay. If semileptonic B decays into the 
four D** states indeed account for less than 2% of the B width, then about another 2% of B decays must be 
semileptonic decays to higher mass excitations or nonresonant channels. 

5. Constrain the validity of models constructed to fit the decays to the ground state final state (most predict, for 
example, larger B D2 than B ^ Di semileptonic rate). 



B. i? — > Xsj photon spectrum 

Comparison of the measured weak radiative B X^j decay rate with theory is an important test of the standard 
model. In contrast to the decay rate itself, the shape of the photon spectrum is not expected to be sensitive to 
new physics, but it can nevertheless provide important information. First of all, studying the photon spectrum is 
important for understanding how precisely the total rate can be predicted in the presence of an experimental cut on 
the photon energy JTot , which is needed for a model independent interpretation of the resulting decay rate. Secondly, 
moments of the photon spectrum may be used to measure the HQET parameters A and Ai p^-p^. 

An important observable is 

dxB (1 - xb) 



N Jl-5 dxB , 

(l-a^s) -1 — , (6) 

xb>1-o / dl 

/ dxB 3 — 
Ji-S axB 

where xb — 2E^/mB is the rescaled photon energy. The parameter S = 1 — 2_E™'"/r7iB corresponds to the experimental 
lower cut on Ej , and it has to satisfy S > Aqcd /"ti b ; otherwise nonperturbative effects are not under control. It is 
straightforward to show that 



(I-xb) = + 1 / dx,{l~Xb) — - S{l-6) — — 

XB>i~s TUB \ niB J Ji^s dxb mB lo dxb 



_ + ..., (7) 

Xb—l—S 



where xi, — 2E^/mi,, and dV/dxi, is the photon spectrum in b quark decay. It has recently been computed away 
from the cndpoint {xb = 1) to order alPo Tq = Gj^ \VtbV^*J'^ rn^ / 32tt'^ is the contribution of the tree 

level matrix element of 07(~ SLC^^F'^'^bR) to the B Xs^ decay rate. All terms but the first one on the right- 
hand-side of Eq. (|^) have perturbative expansions which begin at order as . The ellipses denote contributions of order 
(Aqcd/™b)'^, Q;s(AQCD/»7is)^, and terms not enhanced by /3o, but it does not contain contributions of order 
(Aqcd/™b)^ or additional terms of order (Aqcd /'tis). Terms in the operator product expansion proportional to 
^1,2/^1 enter precisely in the form so that they are absorbed in ms in Eq. (^ ]T^ . 

A determination of A is straightforward using Eq. (^. The left hand side is directly measurable, while the quantities 
entering on the right-hand-side are presented in Ref. The CLEO data in the region E^ > 2.1 GeV yields the 
central values A.^^/^^ ~ 270 MeV and Aq,^ ~ 390 MeV. Uncertainties due to the unknown order Aqj^j-j/mj^ terms in the 
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OPE give the short dashed elUpse in Fig. 2 |T^, whose major axis is roughly perpendicular to those from B Xccv. 
This is why it is important to determine A and Ai from both analyses. The potentially most serious uncertainty is 
from both nonperturbative and perturbative terms that are singular as — > 1 and sum into a shape function [3^ ] 
that modifies the spectrum near the endpoint. For sufficiently large S these effects are not important. They have 
been estimated in Refs. ||l^,|l8| using phenomenological models. Whether these effects are small in a certain range of 
S can be tested experimentally by checking if the extracted value of A is independent of d. This would also improve 
our confidence that the total decay rate in the region xs > 1 — 5 can be predicted model independently. 

There have been other important developments for the B Xg'y decay rate. The next-to-leading order computation 
of the total rate has been completed [p3|-p6[, reducing the theoretical uncertainties significantly. It was realized by 
Voloshin that beyond leading order there are terms in the OPE for the decay rate which are suppressed only by 
^qcd/"^c instead of AgQ^/m^ In fact, there is a series of contributions of the form (AgQ^/m^) (AqcoWb/m^)", 
which gives for n = a calculable correction of ST/T = — (C2/9C7)(A2/m^) ~ 2.5%. It is also known that there are 
uncalculable contributions suppressed by but not by Aqqb /'m-b from photon coupling to light quarks, for which 
there is no OPE While no correction larger than a few percent has been identified, a better understanding of 
the nonperturbative contributions would be desirable once the four-quark operators are included. 



III. \Vub\ FROM THE B X^ev HADRON MASS SPECTRUM 



The traditional method for extracting | Vub \ involves a study of the electron energy spectrum in inclusive semileptonic 
B decay. Electrons with energies in the endpoint region vib/'^, > > {iri^ — m|,)/2mB (in the B rest frame, and 
neglecting the pion mass) must arise from b ^ u transition. Since the size of this region is only about 300 MeV, at the 
present time it is not known how to make a model independent prediction for the spectrum in this region. Another 
possibility for extracting \Vub\ is based on reconstructing the neutrino momentum. The idea is to reconstruct pp and 
infer the invariant mass-squared of the hadronic final state, sh = [pb ^q)'^, where q — +pp- Semileptonic B decays 
satisfying sh < rnjy must come from b ^ u transition [p9|-p| . The first analyses of LEP data utilizing this idea have 
been performed recently [Q. 

Both the invariant mass region, sh < rn?D, and the electron endpoint region, ras/'^. > > {rn^ — m|,)/2rnB, 
receive contributions from hadronic final states with invariant masses between niTr and m/j. However, for the electron 
endpoint region the contribution of states with masses nearer to rriD is strongly suppressed kinematically. In fact, in 
the ISGW model ji^ the electron endpoint region is dominated by the tt and the p with higher mass states making a 
small contribution, and this region includes only of order 10% of the B —> X^ev rate. The situation is very different 
for the low invariant mass region, sh < rnjj. Now all states with invariant masses up to m^j contribute without any 
preferential weighting towards the lowest mass ones. In this case the ISGW model suggests the tt and the p comprise 
only about a quarter of the B semileptonic decays to states with sh < rriD, and only of order 10% of the B —> XuCv 
rate is excluded from this region. Consequently, it is much more likely that the first few terms in the OPE provide 
an accurate description of B semileptonic decay in the region sh < ni^ than in the region E^ > (m^ — m^)/2mB- 
Combining the invariant mass constraint with a modest cut on the electron energy will not destroy this conclusion. 

Let us first consider the contribution of dimension three operators in the OPE to the hadronic mass-squared 
spectrum. This is equivalent to b quark decay, and implies a result for dT/dEo dsQ (where Eq — pf, ■ {pb — q)/mh and 
•So = ipb — <lY the energy and invariant mass of the strongly interacting partons arising from the b quark decay) 
that is straightforward to calculate to order a^/3o Even at this leading order in the OPE there are important 

nonperturbative effects that come from the difference between and ms- The most significant effect comes from 
A, and including only it (i.e., neglecting A1.2), the hadronic invariant mass sh is related to sq and Eq by Jl2| 

SH = So + 2AEq + A^ + .... (8) 

Changing variables from (sq, Eq) to {sh, Eg) and integrating Eq over the range 
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- A < £^0 < 7; {sH - 2KmB + m|), (9) 

gives dV/dsH, where h? < sh < Feynman diagrams with only a u quark in the final state contribute at sq — 0, 
which corresponds to the region < sh < Auib- 

Although dT/dsH is integrable in perturbation theory, it has a double logarithmic singularity at sh — AmB- At 
higher orders in perturbation theory, increasing powers of as\n^[{sH — AmB)/m^] appear in the invariant mass 
spectrum. Therefore, dT / dsn in the vicinity of sh = Ams is hard to predict reliably even in perturbation theory. 
(In the region sh ^ Ams nonperturbative effects are also important.) The behavior of the spectrum near sh = AniB 
becomes less important for observables that average over larger regions of the spectrum, such as dT / dsn integrated 
over s// < A^, with significantly greater than Atub- Fig. 3 shows the fraction of B ^ X^eD events in the region 
SH < A^ as a function of A^ for three different values of A. For a certain value of A, Fig. 3 together with Eq. 
can be used to extract \Vub\ from data, up to the nonperturbative effects discussed next. 

In the low mass region, sh ^ Ams, nonperturbative corrections from higher dimension operators in the OPE are 
very important. Just as in the case of the electron endpoint region in semileptonic B decay or the photon energy 
endpoint region in radiative B decay, the most singular terms can be identified and summed into a shape function. 
These shape functions depend on the same infinite set of matrix elements. Since Amt, ~ 2 GeV^ is not too far from 
mj-) , it is necessary to estimate the influence of the nonperturbative effects on the fraction of B decays with sh < mjj . 
It is difficult to estimate this model independently, but upper bounds can be derived on the fraction of i? — > XuCD 
events with sh > A^ assuming that the shape function is positive In the ACCMM model with reasonable 
parameters, the shape function causes a small (i.e., ~ 4% with A ~ 0.4 GeV, and perturbative QCD corrections 
neglected) fraction of the events to have sh > rnjj ||3^ . This suggests that sensitivity to unknown higher dimension 
operators in the OPE will probably not give rise to a large uncertainty in \Vub\ if it is determined from the hadronic 
invariant mass spectrum in the region sh < "m-^- If experimental resolution forces one to consider a significantly 
smaller region, then the sensitivity to higher dimension operators increases rapidly. 

In summary, to extract \Vub\ from dT/dsn with small theoretical uncertainty, one needs to: 

1. move the experimental cut A as close to mu as possible; 

2. determine A (at order a^) with < 50MeV uncertainty. 

Then a determination of \Vub\ with ~ 10% theoretical uncertainty seems feasible. 
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IV. UPSILON EXPANSION 



The main uncertainties in the theoretical predictions for inclusive B decay rates, e.g., B XuCv, arise from the 
ml dependence on the b quark mass and the bad behavior of the series of perturbative corrections when it is written 
in terms of the pole mass. In fact, only the product of these quantities is unambiguous, but perturbative multi-loop 
calculations are most comfortably done in terms of the pole mass. Of course, one would like to eliminate the quark 
mass altogether from the predictions in favor of a physical observable. Here we present a new method of eliminating 
rrib in terms of the T(IS') meson mass Q (instead of tub and A discussed in Sec. II). 

Let us consider the inclusive B X^ev decay rate [|6|. At the scale = mi,, 



1 _ 2.41^ e - 3.22^/3o - 5.18^P'o - ■ ■ ■ - 



(10) 

'& J 

The variable e = 1 denotes the order in the modified expansion. The complete order calculation was done 
recently and the result is about 90% of the a^Po part. In comparison, the expansion of the T(IS') mass in terms 
of mf, E8| has a different structure. 



2m;, 8 1 TT 
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f3o 



(3^^ + 9£ + 2C(3) + ^ + ^).^ + ... ^ (11) 



where £ = ln[^/ {mbasC p)] and Cp = 4/3. In this expansion we assigned to each term one less power of e than the 
power of as, because as we will sketch below, this is the consistent way of combining Eqs. (|l^) and (11). It is also 
convenient to choose the same renormalization scale fi. The prescription of counting [as(mb)]" in B decay rates as 
order e", and [a,j(mb)]" in mx as order g"^i is called the upsilon expansion. Note that it combines different orders 
in the as perturbation series in Eqs. (|l^) and (pT|). 

The theoretical consistency of the upsilon expansion was shown at large orders for the terms containing the highest 
possible power of /3o, and to order including non-Abelian contributions. An explicit calculation using the Borel 
transform of the static quark potential |^9| shows that the coefficient of the order a""*"^ term in Eq. (|ri|) of the form 
{£" + £"~^ + . . . + 1) exponentiates to give exp{£) = ii/{mbasCF), and corrects the mismatch of the power of as 
between the two series. This is also needed for the cancellation of the renormalon ambiguities in the energy levels 
as given by 2mi, plus the potential and kinetic energies ||50| , |5l]| . The infrared sensitivity of Feynman diagrams can 
be studied by introducing a fictitious infrared cutoff A. The infrared sensitive terms are nonanalytic in A^, such as 
(A2)"/2 or A^" InA^, and arise from the low-momentum part of Feynman diagrams. Diagrams which are more infrared 
sensitive, i.e., have contributions (A^)"/^ or A^" In A^ for small values of n, are expected to have larger nonperturbative 
contributions. Linear infrared sensitivity, i.e., terms of order a/A^, are a signal of Aqcd effects, quadratic sensitivity, 
i.e., terms of order A'^ In A^ are a signal of ^qc£, effects, etc. From Refs. ||5^ and follows that the linear infrared 
sensitivity cancels in the upsilon expansion to order (probably to all orders as well, but the demonstration of this 
appears highly non-trivial). 

Substituting Eq. ( [Ill ) into Eq. (|l^) and collecting terms of a given order in e gives |^ 
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T{B ^ X^eD) 



1927r3 V 2 



1 - O.llSe - 0.035BLMe^ - COOSblmc 



9A2 - Ai 

2(mT/2; 



2 



(12) 



where the BLM |^ subscript indicates that only the corrections proportional to the highest power of /Sq have been 
kept. The complete order term is — 0.041e^ The perturbation series, 1 - 0.115e— 0.035BLMe^ - COOSblme^, is 
far better behaved than the series in Eq. (|o|), 1 — 0.17e — O.ISblme^ — 0.12BLMe'^, or the series expressed in terms of 
the MS mass, 1 + 0.30e -t- 0.19BLMe^ + O.OSblmc"^- The uncertainty in the decay rate using Eq. (^2|) is much smaller 
than that in Eq. (p^), both because the perturbation series is better behaved, and because mx is better known (and 
better defined) than rrib. The relation between \Vub\ and the total semileptonic B Xucu decay rate is Q 

IK.I = (3.06 ± 0.08 ± 0.08) X 10- ( ^iE_L^h^) , (13) 

V 0.001 TB J 
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The first error is obtained by assigning an uncertainty in Eq. ( p^ equal to the value of the term and the second 
is from assuming a 100 MeV uncertainty in Eq. (^l]). The scale dependence of \Vub\ due to varying /i in the range 
mb/2 < fj, < 2mb is less than 1%. The uncertainty in Ai makes a negligible contribution to the total error. Of course, 
it is unlikely that B{B X„eP) will be measured without significant experimental cuts, for example, on the hadronic 
invariant mass (see Sec. Ill), but this method should reduce the uncertainties in such analyses as well. 

The B — > XcCP decay depends on both rrif, and rric. It is convenient to express the decay rate in terms of tox and 
Ai instead of rnj, and rric, using Eq. (hlf) and 



mil — nic = niB — rtiD 



Ai 



Ai 



2mB 2mD 

where ttTb = (Sm^. + mB)/4 = 5.313GeV and m/^ ~ (3m/)* + 77io)/4 = 1.973 GeV. We then find 



, 1 ^ 0.533 X [l - 0.096e - 0.029BLMe^ - (O.28A2 + 0.12Ai)/GeV2] 
1927r'^ V 2 ' 



(14) 



(15) 



where the phase space factor has also been expanded in e. For comparison, the perturbation series in this relation 
when written in terms of the pole mass is 1 — 0.12e — 0.06e^ — . . . [M. Equation (|l5|) implies Es] 



iKbl = (41.6 ± 0.8 ± 0.7 ± 0.5) X 10-^ x 77QED 



B{B X^ev) 1.6 ps 



0.105 



TB 



1/2 



(16) 



where r?QED ~ 1.007 is the electromagnetic radiative correction. The uncertainties come from assuming an error in 
Eq. ( p^ equal to the term, a 0.25 GeV^ error in Ai, and a 100 MeV error in Eq. (pl|), respectively. The second 
uncertainty can be removed by determining Ai, as discussed in Sec. II. Other applications, such as for nonleptonic 
decays, exclusive semileptonic decays and B Xg^ photon spectrum were studied in Refs. p5| , p^ . 

The most important uncertainty in this approach is the size of nonperturbative contributions to tot other than 
those which can be absorbed into the b quark mass. By dimensional analysis the size of this correction is of order 
o'^Aqqq, where a ~ l/{mbas) is the Bohr radius of the T. Quantitative estimates, however, vary in a large range, 
and it is preferable to constrain such effects from data. The upsilon expansion yields parameter free predictions for 
(1 — xb)\xb>i-s defined in Eq. (||). The analog of Eq. 



is K 



{1-xb) 



2m B 



1 + O.Olle + 0.019(e2)BLM - (1 - Xb) 



a:5>(2m£(/m-r)(l-i5) 



(17) 
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For > 2.1 GeV this relation gives 0.111, whereas the central value from the CLEO data is around 0.093. Fig. 4 
shows the prediction for (1 — xb)\xb>i-s as a function of S, both at order e and (e^)BLM- The perturbation expansion 
is very well behaved. In Eq. (^) nonperturbative contributions to mx other than those which can be absorbed into 
the b quark mass have been neglected. If the nonperturbative contribution to T mass, Ax, were known, it could 
be included by replacing mx by mx — Ax- For example. Ax = +100 MeV increases (1 — xb) by 7%, so measuring 
(1 — a;^) with such accuracy will have important implications for the physics of quarkonia as well as for B physics. 



V. CONCLUSIONS 

To conclude, let me emphasize the main points, and indicate what data would be important and useful in my 
opinion to address them: 

• Experimental determination of A and Ai from the semileptonic B — > X^eP lepton energy and hadron mass 
spectra will reduce the theoretical uncertainties in \Vcb\ and \Vub\- 

(Need: Double tagged lepton spectrum with smaller errors.) 

• Photon energy spectrum in i? — s- X^^ gives complimentary information on A and Ai, even in the presence of an 
experimental cut on the photon energy. 

(Need: Spectrum with cut on lowered; even a few hundred MeV can reduce the uncertainties significantly.) 

• To distinguish B — > X^eP from B —f XcCv, cutting on the hadronic invariant mass is theoretically cleaner than 
cutting on the lepton energy. 

(Need: Experimental cut on hadron mass as close to mjj as possible, and a precise determination of A.) 

• The upsilon expansion is equivalent to using a short distance b quark mass, but it eliminates mb altogether 
from the theoretical predictions in favor of mx in a simple and consistent manner. It raises several interesting 
theoretical questions, and has many important applications. 
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